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EARLY  TIME  STRUCTURING  OF  VHANES:  PRELIMINARY  RESULTS 


I.  INTRODUCTION 

An  outstanding  problem  in  the  DNA  HANE  community  is  understanding  the 
generation,  evolution,  and  decay  of  plasma  irregularities  associated  with 
high  altitude  nuclear  bursts.  Although  the  focus  of  this  program  has  been 
on  late  time  structure  (e.g.,  the  determination  of  the  freezing  scale 
size),  an  important  aspect  of  the  problem  is  early  time  structure.  It  is 
known  that  structuring  does  occur  at  very  early  times  in  HANEs  (t  <  sec); 
however,  the  exact  structuring  mechanisms  are  not  well  understood,  nor  is 
the  impact  of  early  time  structure  on  late  time  structure  known.  For 
example,  does  the  early  time  structure  act  as  a  seed  mechanism  for  late 
time  structure?  Are  the  scale  lengths  of  late  time  structure  influenced  by 
the  early  time  structure?  Aside  from  structure  questions,  there  is  also 
the  possibility  that  early  time  structure  can  lead  to  the  transport  of 
debris  across  magnetic  field  lines;  for  example,  weapon  debris  may  reach 
much  higher  altitudes  than  expected. 

Most  of  the  research  on  early  time  stucturing  has  focussed  on  bursts 

9 

which  are  super-Alfvenic,  i.e.,  Starfish  or  Checkmate.  For  the  purposes  of 

9 

discussion  we  will  assume  that  the  transition  altitude  for  super-Alfvenic 

9 

to  sub-Alfvenic  bursts  occurs  at  h  ~  1 000  km.  The  actual  height,  of 

course,  depends  upon  ionospheric  conditions  at  the  time  of  burst  (i.e.,  the 

9 

ambient  Alfven  velocity  V^)  and  the  initial  expansion  velocity  of  the 

debris  (i.e.,  Vjq).  We  refer  to  sub-Alfvenic  bursts  as  VHANEs  (Very 

High  Altitude  Nuclear  Explosions)  since  they  occur  for  h  >  1000  km. 

The  initial  work  on  early  time  structuring  was  done  by  Brecht  and 

Papadopoulos  (1979)  who  suggested  that  a  Rayleigh-Taylor  instability  can 

initiate  cross-field  jetting  of  energetic  ions.  Their  analysis  is  based 

upon  the  usual  MHD  assumptions,  i.e.,  3/9t  «  and  «  R,  Lr:  time 

scales  small  compare  to  the  ion  cyclotron  period,  and  length  scales  large 

compared  to  the  ion  Larmor  radius.  They  argued  that  the  instability  can  be 

driven  by  (1)  laminar  acceleration  of  a  group  of  ions  within  the  debris-air 

shock,  and/or  (2)  the  centripetal  force  associated  with  the  curved  magnetic 

field.  For  parameters  relevant  to  Starfish  they  find  that  the  instability 

has  a  growth  rate  y  -  71  sec-*  (or  growth  time  t  -  14  msec)  and  a 

E 

transverse  wavelength  in  the  range  3-30  km.  However,  their  analysis  is  not 
entirely  applicable  to  VHANEs  (bursts  at  altitudes  h  >  1000  km)  because 

Manuscript  approved  April  26,  1990. 
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their  results  are  based  on  relationships  appropriate  for  super-Alfvenic 
plasma  expansions.  For  example,  they  assume  a  highly  compressed  magnetic 
field  in  the  debris-air  coupling  shell,  and  neglect  the  deceleration  of  the 
debris  shell  caused  by  the  'sveep-up'  of  the  ambient  magnetic  field.  These 

9 

considerations  are  reasonable  for  super-Alfvenic  bursts;  however,  they  are 

9 

not  valid  for  sub-Alfvenic  bursts.  Moreover,  recent  theoretical  work  at 

t 

NRL  on  the  structuring  of  sub-Alfvenic  plasma  expansions  [e.g.,  AM PTE 

magnetotail  release  (Bernhardt  et  al.,  1987)  and  the  NRL  laser  experiment 

(Ripin  et  al.,  1987)]  has  found  that  conventional  MHD  theory  is  not  valid, 

and  that  a  modified  MHD  theory  must  be  used  which  includes  the  Hall  current 

(Hassam  and  Huba,  1987,1988;  Huba  et  al.,  1987,  1989).  Of  course,  hybrid 

codes  (Thomas  and  Brecht,  1988)  and  particle  codes  (Vinske,  1989)  are  also 

appropriate  because  they  contain  the  proper  ion  dynamics. 

Very  recently,  Sperling  (1989)  has  argued  that  the  unmagnetized  ion 

Rayleigh-Taylor  instability  (or  the  large  Larmor  radius  instability  as 

referred  to  in  the  NRL  laser  experiment)  is  not  relevant  to  VHANEs.  A  key 

argument  in  this  paper  is  that  the  conventional  Rayleigh-Taylor  growth  rate 

1/2 

is  less  than  the  ion  cyclotron  frequency  (i.e.,  yq  -  (g/Ln)  <  Bj). 
First,  ve  note  that  this  conclusion  is  based  only  on  a  limited  set  of 
parameters  which  are  poorly  known.  Second,  and  more  important,  ve  have 
shown  in  Huba  et  al.  (1989)  and  in  the  present  paper  that  the  criterion  Yq 

<  fi.  is  overly  stringent,  and  that  a  high  frequency  instability  can  occur 

1  1/4 

for  less  stringent  conditions;  namely,  for  yq  >  (me/mj)  Q^.  Thus,  the 

modified  MHD  theory  of  Hassam  and  Huba  (1988),  and  the  recently  developed 

f 

kinetic  theory  of  sub-Alfvenic  expansions  may  apply  to  VHANEs. 

The  organization  of  the  paper  is  as  follovs.  In  the  next  section  ve 
estimate  the  'effective  gravity'  associated  with  HANE  and  VHANE  plasma 
expansions.  These  results  are  based  upon  simplistic  models  and  provide  a 
good  starting  point  for  more  detailed  calculations  and  numerical 
investigations.  In  Section  III  ve  describe  the  kinetic  theory  of  Huba  et 
al.  (1989)  as  it  pertains  to  the  stability  of  sub-Alfvenic  VHANE 
expansions.  In  Section  IV  ve  apply  the  results  of  Sections  II  and  III  for 
bursts  at  altitudes  h  -  400  km,  1,000  km,  and  10,000  km.  Finally,  in 
Section  V  ve  summarize  our  results. 
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II .  DECELERATION  MODELS 
A.  Mass  Pick-up 

Ve  first  estimate  the  deceleration  of  the  debris-air  coupling  shell 
caused  by  direct  momentum  coupling  to  the  background  plasma.  Ve  estimate 
the  deceleration  of  the  shell  by  assuming  a  momentum  conserving  snowplow 
model.  Ve  assume  that 


^d  -  Vdo  (1> 

where  M  is  the  mass,  is  the  velocity,  and  the  subscript  0  denotes 

initial  values.  Assuming  a  spherical  expansion  we  take 


M 


M0  ♦  j"*3Va 


(2) 


where  R  is  the  radius  of  the  shell,  and  n  and  m  are  the  density  and  mass 

3  a. 

of  the  background  plasma,  respectively.  Ve  substitute  (2)  into  (1)  and 
obtain 


(3) 


1/3 

where  R^  is  the  equal  mass  radius  given  by  R^  -  (3Mo/Annama)  .  Taking 
the  first  derivative  of  (3)  we  obtain  the  effective  gravitational 
acceleration  gM  associated  with  the  deceleration 


*M 


3V  ^  ~ 

JVdO  -3 


*s 


l1 


(4) 


From  (4)  we  find  that  the  maximum  inertial  deceleration  occurs  when  R  * 
0.66Rm  so  that 


(5) 


where  g^m  denotes  the  maximum  'effective'  gravity  (i.e.,  the  deceleration 
dVd/dt). 
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B.  Magnetic  Deceleration 
Ve  estimate  the  'stopping' 


radius  of 


magnetically  confined 


expansion  by  equating  the  initial  kinetic  energy  of  the  plasma  with  the 

3 

magnetic  energy  in  a  volume  (4/3)nRg  where  Rfi  is  defined  as  the  magnetic 
confinement  radius, 


r  Voo 


iHW) 


where  Mq  is  the  mass  of  the  plasma,  V^q  is  the  initial  debris  expansion 
velocity,  and  Bq  is  the  ambient  magnetic  field.  From  (6)  we  obtain 


*B  ■ 


3Hovdo/Bo2) 


Ve  compare  this  distance  with  the  equal  mass  radius  R^.  For  the  expanding 
plasma  to  be  confined  magnetically  we  require  <  RM>  This  leads  to  VdQ 
<  VAa  where  VAa  -  B0/(4nnama)1/^  is  the  Alfven  velocity  in  t.ne  ambient 
plasma.  Vhen  Vd0  «  VAa  the  expanding  plasma  is  stopped  before  it  sweeps 
up  very  much  background  plasma. 

Ve  estimate  gfi,  the  deceleration  caused  magnetic  field  sweep-up, 
using  conservation  of  energy.  We  write 


7  MnV/(t)  +  5~~  F  RJ(t)  .  £  m  V, 


0  dO 


where  the  LHS  of  (8)  is  the  sum  of  the  kinetic  energy  of  the  expanding 
debris  and  the  swept-up  magnetic  energy  at  time  t  and  position  R,  and  the 
RHS  is  the  energy  at  t  >  0.  Ve  solve  (8)  for  V<j(t)  and  obtain 


vd(t)  « 


Ve  take  the  time  derivative  of  (9)  to  obtain 


ft  fad. 

gB(t>  "  “  dt  *  2Mq  R  (t)  "  ZRg  ^ 


where  we  have  made  use  of  the  fact  that  -  dR/dt  and  the  definition  of 
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Rg.  Ve  note  that  the  maximum  deceleration  occurs  at  R  -  Rg,  the  radius  of 
maximum  expansion. 

C.  Curvature  Acceleration 

f 

1.  Super-Alfvenic  Expansions 

2 

We  note  that  the  curvature  force  is  given  by  Fc  *  B  /4rR  so  that 
the  acceleration  gg  is 

gc  *  B2/4nRnTmT  (11) 


vhere  Op  «  nfl  +  n^  and  m^ 
use  the  relationship  B2 
(3),  ve  rewrite  (11)  as 


m  +  m..  In  the  super-Alfvenic  regime  one  can 

Qa  Q 

Vj  4rtnama  (Brecht  and  Papdopoulos,  1979)  and 


V2  R^-2  n  m 

Vd0  .  R  .  nd“d 

gC  *  R  1  +  r3  1  +  n  m 


Ve  take  n^m^  »  nama  (Clark,  private  communication),  which  is  typical  of 
Starfish  for  R  <  Rg,  and  find  that  (12)  becomes 


,  y.  4  [,  *ir2 

*c  ■  "A  *  r2J 


We  will  assume  na  is  roughly  constant,  but  take  nd  to  be  a  decreasing 
function  of  radius  as  follows 

3  M0  1 

- 5 - 5-  .  (14) 

d  8,1  d  AR(3R2  +  AR2) 

Substituting  (14)  into  (13)  ve  find  that 


Vd0  f.  e3V 

R 


2AR(3R2  +  AR2} 


vhere  ve  have  made  use  of  the  definition  of  Rg. 
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2.  Sub-Alfvenic  Expansions 

The  assumptions  upon  which  (13)  is  based  are  not  valid  for 

9 

sub-Alfvenic  expansions.  In  particular,  the  magnetic  field  is  not 

2  2 

compressed  (i.e.,  the  relationship  B  -  4nn&ma  is  no  longer  valid)  so 
that  B  -  Bq  where  Bq  is  the  ambient  magnetic  field.  Substituting  (14) 
into  (11)  we  find  that  the  curvature  acceleration  in  the  sub-Alfvenic 
regime  is  given  by 


w0  r  <3r2  +  **2> 


V  0  0 

dO  AR  3R  +  AR 


We  note  that  gc  reaches  a  minimum  at  R  =  AR//3. 

D.  Quantitative  Results 

2 

We  present  Fig.  1  which  is  a  plot  of  g  (cm/sec  )  vs.  R  (km)  for  a  1 

MT  burst  at  altitudes  h  -  400  km,  1,000  km,  and  10,000  km.  The  parameters 

5  -3 

used  are  the  following:  h  -  400  km  fn  «=  5.0  x  10  cm  and  m  -  16  m 

s  23®  P 

(oxygen  background)];  h  «  1,000  km  {n  *  1.0  x  10  cm  and  m  -  m 

a  o  a  P 

(hydrogen  background)];  and  h  =  10,000  km  ]n  -  1.0  cm  and  m  »  m 

a  a  d 

(hydrogen  background)].  In  all  cases  we  have  taken  VdQ  *  2.0  x  10 

cm/sec,  Mq  «  5.0  x  10^  gm,  m^  =  28  mp  (aluminum  debris),  and  Z  *  2.  These 
parameters  are  shown  in  Table  I  along  with  the  relevant  expansion 
parameters. 

9 

For  h  *  400  km  we  note  that  the  Alfven  Mach  number  is  MA  =  V^O^Aa  * 
9.43  so  that  R^  <  Rg,  as  expected.  Thus,  the  h  *  400  km  burst  is  super- 
Alfvenic;  we  use  (4)  and  (15)  to  calculate  the  mass  and  curvature 
decelerations,  respectively.  The  curvature  acceleration  (g^)  dominates 

*  200  km.  For 

R  >  R^/2  the  two  accelerations  are  comparable.  We  note  that  the  total 
effective  gravity  (gT  «  gc  +  gH)  is  reasonably  large  throughout  the 
expansion  to  RM?  gT  -  10^  cm/sec3. 

For  h  •  1,000  and  10,000  km  we  find  that  the  Alfven  Mach  number  is  M. 

-1  -2  A 
-  1.36  x  10  and  4.71  x  10  ,  respectively.  These  bursts  are  sub- 

9 

Alfvenic  and  therefore  Rg  <  RM.  Thus,  we  use  (10)  and  (16)  in  calculating 

gB  and  gc,  respectively.  The  burst  at  h  «  1,000  km  expands  to  a  radius  Rg 

a  1000  km.  Early  in  the  expansion  phase  we  see  that  g„  >  g_  but  that  the 

6  t  2  “ 

magnitude  of  gc  is  relatively  small,  gc  ~  5  x  10  cm/sec  .  Later  in  the 


the  mass  pick-up  deceleration  (gH)  until  R  *  R^/2  where  R 


6 


expansion  phase  gg  dominates  over  gc  and  becomes  reasonably  large,  gfi  > 

108  cm/sec'.  The  burst  at  h  ■  10,000  km  expands  to  a  radius  Rg  «*  5000  km. 

In  this  situation  ve  note  that  gg  >  gc  throughout  the  expansion,  unlike 

the  h  -  1,000  km  burst.  However,  the  magnitude  of  gn  is  somewhat  smaller 

7  a  2 

than  h  -  1,000  km  burst;  ve  see  that  g„  >  10  cm/sec  for  R  >  2000  km, 

8  2° 

reaching  a  maximum  value  gg  =  10  cm/sec  at  R  -  Rg.  The  reduction  in  the 
value  of  gg  as  h  increases  is  simply  because  the  ambient  magnetic  field 
strength  is  decreasing,  all  other  parameters  being  equal. 


III.  STABILITY  THEORY 

Ve  now  explore  the  stability  of  the  expanding  (and  decelerating) 
debris  plasma.  A  detailed  kinetic  theory  of  the  stability  properties  of 

9 

sub-Alfvenic  plasma  expansions  in  finite  0,  collisional  plasmas  has  been 
presented  in  Huba  et  al.  (1989).  Ve  will  not  reproduce  the  derivation  of 
the  dispersion  equation  here  but  refer  the  interested  reader  to  Huba  et 
al.  (1989)  for  details. 

A.  Equilibrium 

The  plasma  configuration  and  slab  geometry  used  in  the  analysis  are 
shown  in  Fig.  2.  Ve  consider  a  plasma  with  a  density  profile  nQ(x)  such 
that  3nQ/3x  <  0;  we  include  a  magnetic  field  B  *  Bq  ez,  an  ambient 
electric  field  E  -  -Eq  &x,  and  a  gravitational  acceleration  g  *  g  #x.  The 
gravitational  acceleration  can  be  associated  with  the  transverse 

9 

deceleration  of  a  sub-Alfvenic  plasma  expansion  in  a  magnetic  field  (i.e., 

9 

g  *  -dVi0/dt),  a  super-Alfvenic  expansion  in  a  background  plasma,  and/or 
the  curvature  acceleration,  as  described  above.  For  example,  in  Fig.  2, 
one  can  imagine  a  plasma  shell  decelerating  as  it  moves  in  the  +x- 
direction.  Ve  also  assume  that  the  time  scale  of  the  instability  is  much 
slower  than  the  electron  cyclotron  frequency  (i.e.,  o>  «  Cg).  Ve  allow 
the  time  scale  of  the  instability  to  be  arbitrary  in  relation  to  the  ion 
cyclotron  period. 

The  electron  momentum  equation  for  this  situation  is  given  by 


0 


+  n 
c~e0 


x 


(17) 


From  (17)  it  is  easily  shown  that  the  equilibrium  electron  drift  is  given 
by 
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2. 0  -  <VE  +  Vde>  ey 


where  Vg  «  cEq/Bq  is  t*ie  1  x  1  drift  velocity  and  *  cTg/eBQLn  is 
electron  diamagnetic  drift  velocity  where  Lr  -  |  Bln  nQ/Bx|“*. 

The  ion  momentum  equation  is  given  by 


Tj  7n, 


n  eZ  ,,,  1„  n  .  i  ¥,'0 

v~0  c~iO  ~0  m^  n^ 


From  (19)  the  equilibrium  ion  drift  is  found  to  be 

V  ,  (ffo  _  £»  _  L.) ; 

~i0  2.  eZBn  L  J  y  ' 


The  components  of  VQ  are  the  E  x  B  drift  [Vg  «  cEq/Bq],  the  gravitational 
drift  [Vg  -  g/fijlt  and  the  ion  diamagnetic  drift  [Vdi  -  cTj/eZBgl^J. 

In  order  to  satisfy  the  equilibrium  condition  V  •  nyiQ  =  0  we  choose 
to  work  in  the  ion  rest  frame,  i.e.,  yiQ  *  0.  Thus,  we  take 

cE0  g  cTi  1 

B  *  q  +  eZB  L  (21) 

B0  8i  eZB0  n 

The  ions  are  electrostatically  confined  in  the  unmagnetized  limit. 

B.  Dispersion  Equation 

For  the  parameters  given  in  Table  I  we  note  3  (=  8raiT/B  )  «  1,  ven 
«  Qe,  and  \>in  «  Qj.  Thus,  we  will  use  the  low  3  (i.e.,  electrostatic), 
collisionless  dispersion  equation  derived  in  Huba  et  al.  (1989).  In  this 
limit  the  dispersion  equation  is 


D(a>,ky) 


W _  2 «  .  2w 

1  +  +  -^2  (1  +  iwGj)  + 


W1  ~  h>2r0(be) 


vhere 
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OB  iC  V 

Gd  -  dt  exp[i(u  +  kyVdi)t  +  §kJp2(cos  8^-1)  +  t-^^sin 


«  -  kyV£,  »2 


“  '  kyVE  "  kyVde’  "pi2  *  and  upe2 

tr  t  r\  tr  /-»tot  l  i-  ^  .  4-  /  *"» 


4nneVme,  V£  -  Vg  +  Vdi,  Vg  -  g/Q^  Vdi  .  cT^eZBL,,  bg  -  ky‘p%‘/2,  p  - 

kyPe’  pe  ■  Ve/Se’  Va2  *  2Ta/ma»  ro(x)  =  exP<-x>Io<x> »  and  Io(x>  is  the 
modified  Bessel  function  of  order  0.  This  dispersion  equation  is  valid 

for  (o  «  Oj,  <*>  -  Sj,  and  w  »  B^;  thus,  it  describes  the  Rayleigh-Taylor 
instability  in  both  the  magnetized  and  unmagnetized  ion  limits. 

1.  Figh  Frequency  Limit  (w  »  8^) 

We  first  simplify  (22)  by  considering  the  high  frequency  (y  >  Qj 
where  y  is  the  growth  rate),  short  wavelength  (kypj  »  1),  cold  electron 
limit  (Tfi  •+  0).  In  this  limit  we  note  that  Gj  *  -i^Z(^)/w  where  ^  * 
w/kyvit  and  Ig(be)  =  1  -  bg.  The  dispersion  equation  then  becomes 


D(o>fky) 


2<.  2w  .  k  V.. 

1  +  +  1  ^ 2  (1  +  5jZ(^))  +  ~2^2  (lo'-  ItV 

k«vi  k,,v4  y  E 


We  further  simplify  (23)  by  assuming  the  ions  to  be  either  cold 

(«/k v.  »  1)  or  warm  («/k  v.  «  1),  which  corresponds  to  the  strong  drift 
y  i  y  1 

velocity  regime  (V£  »  vp  or  the  weak  drift  velocity  regime  (V£  «  v^, 
respectively. 

a.  Cold  Ion  Limit  («/kyv^  »  1) 

We  first  consider  the  cold  ion  limit:  w/kyV.^  >>  1  {the 
strong  drift  regime  (i.e.,  V£  »  v^)J.  We  expand  the  plasma  dispersion 

function  in  (23)  in  its  asymptotic  limit.  The  dispersion  equation  then 
becomes 

2  2  2 

2®„l  kvVrf< 

D(w,k  )  *  1  -  —j—  +  ^  u  v  +  *  0  (24) 

7  «  kyV‘  "  "  kyVE  82 

We  can  rewrite  (24)  in  the  following  form 


<“2  -  “lh><“  -  kyVE>  * 


2  2 
"lh" 

SikyLn 
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2  2  1/2 

where  -  0^/(1  +  w ^  /Q^  )  is  the  lover  hybrid  frequency.  Equation 

(25)  is  very  similar  to  the  dispersion  equation  derived  by  Krall  and 

Liewer  (1971)  for  the  lower-hybrid-drift  instability  (see  Eq.  (7)  of  their 

paper);  in  the  limit  g  -»  0  it  agrees  exactly  with  their  dispersion 

equation.  Following  the  analysis  of  Krall  and  Liewer  (1971)  one  can  show 

that  the  maximum  growth  rate  of  the  instability  is  rM  *  with  »  kyVg 

a  From  (25)  it  is  clear  that  the  instability  is  fluid-like  and  is 

caused  by  the  coupling  of  a  lower  hybrid  wave  (co  =  w^)  with  a  drift  wave 

(u  :  k  V.)  in  an  inhomogeneous  plasma  (i.e.,  finite  L  ). 
y  t  n 

Equation  (25)  yields  a  relatively  simple  solution  in  the  low 

frequency  limit,  i.e.,  w  <<  It  is  given  by 

«2  -  Q.kyLjjW  +  k2gLn  =  0  (26) 

where  we  have  assumed  >>  V^.  Equation  (26)  has  previously  been 

derived  from  fluid  theory  by  Hassam  and  Huba  (1987,  1988).  The 

eigenfrequency  is  given  by 


“  ■  I  VyL„  *  I  «aiV»)2  -  4gLnk5|1/2  '  <27> 

2 

Instability  occurs  for  g/L  >  Q.  /A.  In  the  limit  of  large  g  one  finds 

1/2  n  1 
that  Y  =  ky(gLn)i/z. 

b.  Varm  Ion  Limit  (w/kyVj  <<  1) 

Ve  now  consider  the  varm  ion  limit:  «/kyv^  «  1  [the  weak 
drift  regime  (i.e.,  Vg  «  Vj)J.  For  this  situation  it  is  found  that  the 
instability  is  kinetic  in  nature;  the  unstable  waves  grow  because  of 
inverse  Landau  damping  of  the  ions.  Ve  expand  the  plasma  dispersion 
function  in  (23)  in  the  small  argument  limit.  The  dispersion  equation 
then  becomes 
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(29) 


j2  1.2  1.  y  |2 

JCM  km  K  vdl  ^  1/jlw 

°-1%7*k7^vi^vr 


where 


k2  rii(1,5sri  . 

°e 

Assuming  that  «r  »  y,  the  real  frequency  is  determined  from  D(wr,ky)  *  0 
and  we  find  that 


«  =  k  V, .  — s— — s-  +  k  V  . 

r  y  41k*  ♦  ^  y  g 


From  (30)  we  see  that 


kV, 


to  s  —  +  k  V 

r  kj  yg 


for  ky  «  kM 


(31a) 


•r  =  VVdi  +  Vg> 


for  ky  »  kM 


(31b) 


There  is  a  significant  difference  in  the  ky  dependence  of  the  real 

frequency  betveen  the  ion  diamagnetic  drift  and  the  gravitational  drift  in 

the  long  wavelength  regime,  k  «  1c..  For  the  case  of  the  lower-hybrid- 

y  n  _ 

drift  instability  (V^^  »  V^)  it  is  found  that  »r  *  ky  ,  while  for  the 
unmagnetized  ion  Rayleigh-Taylor  instability  (V  »  Vdi)  o>r  «  ky. 

The  growth  rate  is  de*  rmined  by  the  expression  y  * 
-Di(»r*ky)/(3Dr/3»r)  yhere  the  subscripts  r  and  i  refer  to  the  real  and 

imaginary  components  of  D,  respectively.  We  find  that 


(  kM  \2tVM  kv  V 

(r2  rr)  (“  .2  1  .2  *  vf. 
\  *  Si  1  ky*  Si  1 


)  JT 


For  the  case  V  ■  0,  (30)  and  (32)  reduce  to  the  expressions  for  co 

O  ^ 

and  r  derived  in  Davidson  et  al.  (1977)  for  the  lower-hybrid-drift 
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instability.  It  can  be  shovn  that  the  growth  rate  maximizes  at  ky  -  k^j 
the  grovth  rate  and  real  frequency  at  k^  -  k^  are  given  by 


33 

8  2  “lh 

vi 


%  *  I  Vdi 


(34) 


In  the  opposite  limit,  when  the  gravitational  drift  dominates  (i.e., 
Vg  »  Vjj)  we  find  that  the  maximum  growth  rate  occurs  at  ky  •  k^Z/3;  the 
corresponding  values  of  the  grovth  rate  and  real  frequency  are  given  by 


One  interesting  difference  between  the  grovth  rate  in  the  strong  drift 

limit  and  weak  drift  limit  is  the  dependence  on  the  density  gradient  scale 

1/2 

length.  In  the  former  case  we  found  that  y  *  L  (see  (27)),  while  in 

-1  n 

the  latter  case  y  *  L„  (see  (34)).  The  requirement  for  the  unmagnetized 
ion  Rayleigh-Taylor  instability  to  be  unstable  in  the  kinetic  regime  is  y 

>  8. ;  from  (34)  we  find  that  the  turn-on  condition  is  roughly  g/L  > 

1  1/2  2  n 
(me/mj)  8j  .  This  criterion  is  considerably  less  stringent  than  that 

for  the  unmagnetized  ion  Rayleigh-Taylor  instability  in  the  fluid  limit 

(i.e.,  g/Ln  >  8i2/4). 

2.  Low  Frequency  Limit  («  <<  8^) 

We  now  consider  the  low  frequency  («  «  8^),  short  wavelength 
(kyPj  «  1)  limit.  In  this  limit  we  note  that 


L)d  -  bj) 


and  the  dispersion  equation  becomes 
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0 


(37) 


where  b1  -  V.^o^/2.  In  the  limit  that  «  j2/*^2  »  Wpe2/Qe2  >>  1  ll  is 
easily  shown  that  (37)  reduces  to 


<"2  -  (fcyVg  ♦  kyVdi)“  *  fr  -  0  (38) 

which  is  the  usual  dispersion  equation  for  the  Rayleigh-Taylor  instability 

in  the  limit  of  magnetized  ions.  In  the  limit  «  »  k  V  ,  k  V ,  the 

1/2  y  8  y  ai 

growth  rate  is  simply  given  by  y  *  (g/Ln)  '  The  instability  is 
stabilized  when  (kyV^  +  k^V^)2  *  *g/Ln;  *s  tlle  so-calle(*  finite 

Larmor  radius  stabilization  of  the  Rayleigh-Taylor  instability  (Roberts 
and  Taylor,  1962). 


IV.  APPLICATION  TO  VHANES 

Ve  now  apply  the  results  of  the  stability  analysis  presented  in 
Section  III  to  the  nuclear  burst  parameters  presented  in  Section  II.  Ve 
first  write  down  the  turn-on  conditions  for  the  instabilities  discussed  in 
the  previous  section.  Ve  cast  these  conditions  in  terms  of  the  growth 
rate  of  the  conventional  Rayleigh-Taylor  instability,  i.e.,  yn  * 

1  /O  ^ 


Instability  Turn-on  Criterion 

Unmagnetized  Ion  Rayleigh-  yn  >  B./2 

Taylor:  Fluid  Limit 

1/4 

Unmagnetized  Ion  Rayleigh-  yn  >  (m  /m. )  8. 

Taylor:  Kinetic  Limit 


Magnetized  Ion  Rayleigh-  y_  >  (k  V  +  k  V,,)/2 

Taylor:  MHD  Limit  u  y  g  y 


In  Fig.  3  we  plot  y^/B^  vs.  R  for  the  burst  parameters  shown  in  Fig. 
1.  In  calculating  y^  we  have  taken  the  density  gradient  scale  length  to 
be  L  «  50  km.  Ve  also  indicate  the  turn-on  values  for  the  unmagnetized 
ion  Rayleigh-Taylor  instabilities:  y^/fi^  >  0.5  the  fluid  mode  is  unstable, 
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and  Yq/S^  >  0.066  Che  kinetic  mode  is  unstable.  The  magnetized  ion 

Rayleigh-Taylor  instability  will  be  unstable  for  sufficiently  small 
vavenumbers  (i.e.,  k^p^  «  1).  For  the  burst  at  altitude  h  -  400  km 
(Starfish),  we  note  that  -  °i0^Vd0/VAa^*  +  r^rm^)~*  which  is  based  on 
the  relationship  -  4nnamaV^ .  Thus,  we  account  for  the  compression  of 
the  magnetic  field  in  the  debris  shell  which  occurs  in  a  super-Alfvenic 
expansion.  Ve  see  that  Yq/Q^  <  0.02  for  R  <  so  that  the  debris  shell 
is  stable  to  the  unmagnetized  ion  Rayleigh-Taylor  instabilities.  For  the 
burst  at  h  >  1,000  ka  we  find  that  the  kinetic  instability  may  go  unstable 
at  R  *  R0.  However,  the  burst  at  h  =  10,000  km  is  strongly  unstable  to 
both  the  kinetic  instability  (which  turns-on  at  R  «  400  km)  and  the  fluid 
instability  (which  turns-on  at  R  =  4000  km). 

An  important  consideration  in  the  above  discussion  is  the  sensitivity 
of  the  stability  criterion  to  the  density  gradient  scale  length.  For 

example,  if  we  had  chosen  Ln  »  5  ku  for  the  burst  at  h  *  400  km,  then  the 

curve  Yq/Q^  would  be  increased  by  a  factor  of  3  and  the  kinetic 
instability  would  be  unstable  near  R  a  R^.  On  the  other  hand,  if  the 

density  gradient  scale  length  were  Lfi  *  100  km  for  the  burst  at  h  >  1,000 

km,  then  the  expansion  would  be  predicted  to  be  stable  to  the  unmagnetized 
ion  Rayleigh-Taylor  instability.  Moreover,  the  density  gradient  scale 
length  will  probably  not  be  constant  during  the  expansion;  clearly, 
computer  simulations  are  needed  to  better  determine  the  thickness  of  the 
debris  shell  as  it  expands  and  decelerates. 

Finally,  we  present  Fig.  4  which  is  a  plot  of  y/Bj  versus  kypj.  We 
solve  (22)  numerically  for  a  set  of  parameters  relevant  to  the  burst  at  h 
■  10,000  km  at  an  expansion  radius  R  =  5,000  km.  Ve  consider  the 
following  parameters:  V^/v^  *  4.0,  V^/v^  *  0.014,  Ln  =  50  km,  Te/Tj  * 

1.0,  T.  *  T,  «  T  *  50  ev,  p.  -  0.0,  co  /B  «  10.0,  Z  =  2,  and  g,  ■  28 

l  a  a  i  pe  e  o 

(aluminum).  The  two  branches  of  unstable  modes  are  apparent.  In  the  low 
frequency  (y  <  S^),  long  wavelength  limit  (kyPj  «  1)  we  obtain  the 
conventional  magnetized  ion  Rayleigh-Taylor  instability;  for  kyp^  «  1  we 

find  that  the  growth  rate  asymptotes  to  the  standard  growth  rate  yn  « 

1/2  u 
(g/Ln)  «  0.335  Qj.  The  local  approximation  breaks  down  for  kyLn  1  1 

which  occurs  for  k^p^  a  0.028  so  that  the  results  presented  are  valid.  As 

kyp^  increases  above  k^p^  -  0.15  we  see  that  the  mode  becomes  stable 

because  of  finite  Larmor  radius  effects  and  the  real  frequency  associated 
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with  the  gravitational  drift  vave.  In  the  high  frequency  (y  >  P, ) ,  short 
vavelength  limit  (kyPj  »  1)  ve  find  the  unmagnetized  ion  Rayleigh-Taylor 
instability  vhich  is  driven  by  a  coupling  betveen  the  gravitational  drift 

vave  and  the  lover  hybrid  vave.  The  grovth  rate  maximizes  for  b  e  b,, 

1/2  r  in 
vhere  »  ( 2eQ^ )  is  the  lover  hybrid  frequency.  For  the  parameters 

chosen,  ve  find  that  the  grovth  rate  maximizes  at  k^p^  a  40  vhich 

corresponds  to  a  vavelength  X  a  200  m.  On  the  other  hand,  if  ve  were  to 

assume  that  the  dominant  vavelength  is  X  =  Ln,  this  vould  correspond  to 

kyp^  a  0.176  vhich  is  stable  for  the  parameters  chosen. 

Ve  add  that  ve  have  chosen  a  single  set  of  parameters;  based  on  these 
parameters  our  results  suggest  that  the  unmagnetized  ion  Rayleigh-Taylor 
instability  (or  the  large  Larmor  radius  instability)  can  play  a  dominant 
role  in  the  structuring  of  the  debris  shell  for  a  VHANE.  Hovever,  the 
results  are  sensitive  to  several  parameters,  in  particular  the  density 
gradient  scale  length,  and  simulation  studies  are  needed  to  better 
determine  the  macroscopic  parameters  associated  vith  a  VHANE  (Thomas  and 
Brecht,  1988;  Vinske,  1989). 

V.  DISCUSSION 

The  stability  of  the  debris  shell  of  a  nuclear  burst  at  altitudes  h  > 
400  km  is  investigated.  A  set  of  relatively  simple  expressions  are 
derived  to  estimate  the  effective  gravitational  acceleration  associated 
vith  the  deceleration  of  the  shell  (by  mass  'pick-up'  and  magnetic  field 
'sveep-up')  and  the  curvature  of  the  magnetic  field.  Ve  then  present  a 
stability  analysis  based  on  the  recent  kinetic  theory  developed  by  Huba  et 
al.  (1989).  In  particular,  the  'turn-on'  conditions  for  the  unmagnetized 
ion  Rayleigh-Taylor  instability  are  derived  for  both  the  fluid  and  kinetic 
regimes,  as  veil  as  the  finite  Larmor  radius  stabilization  criterion  for 
the  magnetized  ion  Rayleigh-Taylor  instability.  Ve  apply  these  results  to 
1  NT  bursts  at  altitudes  h  -  400  km,  1,000  km,  and  10,000  km.  Ve  find  the 
burst  at  h  ■  400  km  to  be  stable  to  the  unmagnetized  ion  Rayleigh-Taylor 
instability;  the  burst  at  h  -  1,000  is  marginally  unstable  to  the  kinetic 

Instability;  and  the  burst  at  h  -  10,000  km  is  strongly  unstable  to  both 

the  kinetic  and  fluid  unmagnetized  ion  Rayleigh-Taylor  instabilities.  On 
the  other  hand,  the  magnetized  ion  Rayleigh-Taylor  instability  is  unstable 
at  all  altitudes  and  may  cause  structuring  of  the  debris  shell  depending 
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on  the  growth  rate  of  the  instability.  A  critical  parameter  in 
determining  the  stability  properties  of  the  debris  shell  is  the  density 
gradient  scale  length  (or  shell  thickness).  It  is  recommended  that 
detailed  2D  and/or  3D  computer  simulations  be  performed  to  quantify  the 
macroscopic  parameters  associated  with  the  expansion  of  the  debris  shell 
(Thomas  and  Brecht,  1988;  Vinske,  1989) 
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TABLE  Is  BURST  PARAMETERS 


Initial  expansion  velocity  (VdQ) 
Mass  of  burst  (MQ) 

Ambient  magnetic  field  (G) 

Debris  atomic  number  (ud> 
Background  atomic  number  (ya) 
Ambient  density  (n&) 

Debris  temperature  (T^) 

Debris  thermal  velocity  (vd) 

Debris  charge  state  (Z) 

Density  gradient  scale  length  (Lr) 
Ambient  cyclotron  frequency  (S^ 
Directed  Larmor  radius  (V^/8^) 
Ambient  Alfven  velocity  (Vda) 
Alfven  Mach  number  (M^) 

Equal  mass  radius  (R^) 

Magnetic  confinement  radius  (Rg) 


Burst 

Altitude  h  (km) 

400 

1,000 

2000  km/ sec 

2000  km/sec 

5.0  x  105  gm 

5.0  x  103  gm 

0.28 

0.21 

28 

28 

16 

1 

5.0  x  105  cm-3 

1.0  x  103  cm"3 

50  ev 

50  ev 

18.5  km/sec 

18.5  km/sec 

2 

2 

50  km 

50  km 

188  rad/sec 

146  rad/sec 

10.6  km 

13.7  km 

2  x  102  km/sec 

2  x  104  km/sec 

9.43 

0,136 

2.0  x  102  km 

4.1  x  103  km 

9.1  x  102  km 

1.1  x  103  km 

10,000 
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1E7  -U4 - - - * - - - * - • — ■  ■>  — 1 - 

10  30  100  200 

R  (km) 

Fig.  1  —  Plot  of  g  (cm/ sec2)  vs  R  (km)  for  a  1  MT  burst  at  altitudes  h  =  400  km,  1,000  km,  and  10,000  km.  The  parame¬ 
ters  used  are  the  following:  h  =  400  km  [n,  =  5.0  x  10s  cm'3  and  m,  =  16  m^  (oxygen  background)];  h  =  1,000  km 
(n,  *  1.0  x  103  cm'3  and  m,  =  (hydrogen  background)];  and  h  =  100,000  km  [n,  =  1.0  cm_3  and  m,  =  m,,  (hydro¬ 
gen  background)].  In  all  cases  we  have  taken  =  2.0  x  10*  cm/sec,  Mo  =  5.0  x  103  gm,  rr^  =  28  nip  (aluminum 
debris),  and  Z  -  2.  These  parameters  are  shown  in  Table  I  along  with  the  relevant  expansion  parameters. 
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Fig.  3  —  Plot  of  7o/0,  vs  R  for  the  burst  parameters  shown  in  Fig.  1.  In  calculating  -y0  we  have  taken  the  density  gradient 
scale  length  to  be  L,  =  50  km.  We  also  indicate  the  turn-on  values  for  the  unmagnetized  ion  Rayieigh-Taylor  instabilities; 
y0/Q>  >  0.5  the  fluid  mode  is  unstable,  and  y0/Q>  >  0.066  the  kinetic  mode  is  unstable.  The  magnetized  ion  Rayieigh- 
Taylor  instability  will  be  unstable  for  sufficiently  small  wavenumbers  (i.e. ,  k^p,  «  1). 
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I  Magnetized  ion  RT 


I  I  H  tl 

70.0 

Fig.  4  —  Plot  of  y/(l ,  vs  kyp,  for  a  set  of  parameters  relevant  to  the  burst  at  h  =  10.000  km  at  an  expansion  radius  R  = 
5,000  km.  We  consider  the  following  parameters:  Vg/v,  =  4.0,  Vdj/v,  =  0.014,  L„  =  50  km,  Tt/T,  =  1.0, 
T,  =  Td  =  T,  =  50  ev,  0,  =  0.0,  =  10.0,  Z  =  2,  and  ^  =  28  (aluminum). 


/  i 

Unmagnetized  ion  RT 
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